Inverse statistics of smooth signals: the case of two dimensional turbulence 
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The problem of inverse statistics (statistics of distances for which the signal fluctuations are 
larger than a certain threshold) in differentiable signals with power law spectrum, E(k) ~ k~ a , 
3 < a < 5, is discussed. We show that for these signals, with random phases, exit-distance moments 
follow a bi-fractal distribution. We also investigate two dimensional turbulent flows in the direct 
cascade regime, which display a more complex behavior. We give numerical evidences that the 
inverse statistics of 2d turbulent flows is described by a multi-fractal probability distribution, i.e. 
the statistics of laminar events is not simply captured by the exponent a characterizing the spectrum. 
PACS number(s) : 47.27.Eq, 05.45.-a 



Many phenomena in natural science posses complex be- 
haviors over a wide range of spatial and temporal scales. 
Complexity is quantified by the non-Gaussian properties 
of probability distribution functions (pdf) of signal in- 
crements over a given range of scales. Whenever pdf's at 
different scales cannot be superposed by a simple rescal- 
ing procedure one speaks about intermittency Such 
non-trivial rescaling properties may be exhibited by pdf's 
tails or peaks, or both ||. Strongly intermittent and 
rough signals, like those encountered in three dimen- 
sional turbulent flows, are the typical examples of sys- 
tems with non-trivial, say multiaffine, scaling of pdf's 
tails. Other important natural phenomena develop sim- 
ple pdf's tails but non-trivial pdf's peaks. This is the 
case of two dimensional turbulence as it will be shown in 
this letter. Indeed laminar fluctuations, corresponding 
to the events described by the peak of the probability 
distribution, posses non-trivial scaling properties. Re- 
cently, it has been shown that laminar fluctuations of 
rough and multiaffine fields are optimally characterized 
in terms of their exit-distance statistics, also known as 
inverse-statistics 

The aim of this letter is twofold. First we want to 
extend the application of inverse statistics ^|JJ] to the 
case of smooth signals with a given power spectrum, 
E{k) ~ k~ a . In particular, we discuss signals only 
one time differentiable, i.e. 3 < a < 5. For such 
signals, direct statistics, i.e. moments of signal incre- 
ments over a given scale, do not bring any information: 
they are always dominated by the differentiable event, 
v(x + r) — v(x) ~ r, which are always present when the 
spectrum has a slope a > 3. On the contrary, we will 
show that the exit-distance statistics is given by a bi- 
fractal distribution if only one kind of more than smooth 
fluctuation exists. With more than smooth fluctuations 
we mean events where the signal has a local scaling as 
v(x + r) — v{x) ~ r h with h > 1. 

Second, we will apply the inverse statistics analysis to 
the case of two dimensional turbulence. Two-dimensional 
turbulence in the direct enstrophy cascade regime is of 
obvious importance both theoretically and practically to 
understand a variety of different natural processes, e.g. in 



geophysics and astrophysics By applying the exit- 

distance analysis to a set of 2d numerical simulations, 
we will show in a quantitative way that laminar events 
of two-dimensional turbulence posses highly non-trivial 
rescaling properties, revealing a rich (multifractal) struc- 
ture of laminar fluctuations. We also discuss the impor- 
tance of large-scale structures in determining the inverse 
statistics of two dimensional turbulent flows. 

Let us begin by considering a one dimensional signal 
built by fixing its spectrum as E(k) ~ k~ a on all availa- 
ble wave-numbers : 
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with |w(A:)| 2 ~ fc~ Q and 9k random phases, uniformly dis- 
tributed in [0,27r]. If 3 < a < 5 the signal is everywhere 
one-time differentiable : it follows that moments of its dif- 
ference over any increment r always have a differentiable 
scaling, namely 

S p (r) = ({v(x + r) - v(x))P) ~ Cp rP , 

while moments with p < —1 do not exist. In order 
to highlight the role of non-trivial more than smooth 
stochastic fluctuations associated to the spectrum slope 
a, one needs observable sensitive to laminar fluctuations 
like moments of inverse statistics. With inverse statistics 
we mean moments of increments, r(Sv), necessary to ob- 
serve in the signal a forward (backward) exit through a 
barrier Sv. In particular, we fix the height of the barrier, 
5v, and we pick at random a reference point xq. Then, 
we measure the first forward or backward exit, r(Sv), i.e. 
the first point such that \v(xq ± r) — u(a?o)| > & v an d 
we repeat the observations for many xq and for different 
barrier heights. This allows us to define the probability 
distribution for exit events, V(r(Sv)). 
Let us define, the exit moments through a barrier Sv as: 



T (p \Sv) = (r p (5v)), 



(2) 



where the average is taken with respect to the random 
choice of xq Q . Positive moments of exit events prefer- 
entially weight smooth, laminar, fluctuations. Exit-event 
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moments are also called inverse statistics moments be- 
cause of the difference with direct statistics where one 
measure signal increments over a given scale. For the 
prototype smooth signal ([!]), a rigorous estimate of the 
scaling exponents of inverse statistics moments can be de- 
rived as follows. From ([!]), when the spectrum exponent 
is in the window of first order differentiability, 3 < a < 5, 
we may estimate the typical fluctuations as 



v(x + r) - v(x ) ~ d Xo v(x )r + c(x )r h 



(3) 



where we have kept only the two most important scal- 
ing behaviors: 0(r) because of the differentiability and 
0(r h ) from the spectrum exponent. In (^), the exponent 
1 < h < 2 is connected to the spectrum slope by the 
dimensional relation a = 2h + 1, while the function c(x) 
is a continuous function of x. By studying the exit event, 
in the limit of small barrier height, we select with proba- 
bility one the differentiable scaling r(Sv) ~ Sv except for 
those xqS where the first derivative, d x v(xo), vanishes. 
In the latter case, the 0(r h ) term dominates the scal- 
ing behavior. With 3 < a < 5 the first derivative is a 
sclf-affine signal with Holder exponent £ = h — 1, i.e. 
d x v(x + r) — d x v{x) ~ r^, that is it vanishes on a fractal 
set of dimension D = 1 — £ = 2 — h. Therefore, the prob- 
ability to see the sub-dominant term 0(r h ) dominating 
the exit events in (||) is given by the probability to pick 
a point at random on a fractal set with dimension D, i.e. 



P(r ~ (Sv) 1 / 11 ) ~ r y - u = (Sv) 



l-D 



a-i/h 



(4) 



Taking into account both situations, we end with the 
following bi-fractal prediction for inverse statistics mo- 
ments : 



T ip \Sv) ~ Sv Xbt < p \ xm(p) 
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(5) 



From the previous bi-fractal formula, one sees that lam- 
inar, differentiable, fluctuations influence the inverse 
statistics only up to moments of order p = 1 ; for larger 
p, the pdf is dominated by the sub-dominant behavior, 
v(x + r) — v{x) ~ r . In other words, the extrema of 
the signal play the role of singularities in inverse statis- 
tics : close to the extrema, events with much longer exits 
through barriers of order Sv are observed when Sv — > 0. 
In Fig. 1, we numerically check this prediction on a 
one-dimensional signal with a = 4, i.e. with h = 1.5. 
The derivative of such a signal has an Holder exponent 
£ = 0.5, i.e. it is a stationary Brownian motion. As 
shown, the prediction is verified with high accuracy. 

Switching to real signals, we consider two dimensional 
incompressible turbulent flows. As it is well known, 2d 
turbulence is characterized by a forward enstrophy cas- 
cade from forced scales to dissipative ones Q . In the iner- 
tial range, arguments d la Kolmogorov give for the veloc- 
ity spectrum the prediction E(k) ~ fc -3 plus logarithmic 



corrections ||, which is experimentally ]l0| and numeri- 
cally [Tl| ] observed. However, it is also measured a more 
than smooth spectrum E{k) ~ kr a with a > 3 depend- 
ing on the characteristics of the forcing and large-scale 
dissipation p2 13 1. The dependency of inertial range 



statistics from large-scale effects (universality issue) is 
still an open problem in 2d turbulence in the direct cas- 
cade regime. 




FIG. 1. Scaling exponents x(p) f° r the 1-d signal ([!]) with 
a — 4. The dashed line refers to the linear differentiable be- 
havior for p < 1, Xbf(p) — P- The continuous line gives the 
more than smooth behavior, Xbf(f) = {p/h + 1 — 1/h). Mo- 



have been computed using 10 realizations of the 
with 2 17 modes; for each realization 2 12 starting 
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points, xo, have been taken at random. 

The most important feature of the vorticity cascade is 
the presence of a strong interaction between eddies of 
very different scales : such a non-locality (in Fourier 
space) should play a fundamental role in shaping the en- 
ergy spectrum, in particular the strong dependency of 
the spectrum slope from large-scale effects. In order to 
understand the importance of large-scale statistics one 
needs to go beyond the spectrum slope. 

We analyze exit-distance statistics in a series of direct 
numerical simulations of the two-dimensional Navier- 
Stokes equation: 



d t uj + J{oj, ip) = vk q uj - (3 p A- p uj + F, 



(0) 



where lo is the vorticity, tp the stream function, and J the 
Jacobian. We use a standard dealiased pseudo-spectral 
algorithm with periodic boundary conditions, at resolu- 
tions 512 2 and 1024 2 . The large-scale forcing F is Gaus- 
sian, white-in-time, and nonzero only at some character- 
istic wave-numbers kf between 4 and 6. Enstrophy is 
dissipated at small scales with an hyper- viscosity q = 4. 
Energy is removed at large scales to avoid piling up on 
the smallest mode. We performed two sets of numerical 
simulations with different IR draining, with an inverse 
Laplacian with p — 2 case (A) and with p = case (B), 
without observing big differences in the spectrum slope. 
In Fig. 2 (inset) we show the compensated average spec- 
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trum that we observe for case (A) : notice the strong 
influence of coherent structures for low wave-numbers. 
The best fit to the spectrum slope for this case gives 
a = 3.26 ± 0.06, i.e. a more than smooth exponent 
h = 1.13, while for the case (B) we found a — 3.24±0.06. 
The relevant test we wish to perform on the 2d flow 
is comparing the inverse statistics measured on several 
snapshots of the direct numerical simulations with the 
inverse statistics obtained after randomization of all ve- 
locity phases on the same frames. We measured the mo- 
ments of exit events using both transversal and longitu- 
dinal velocity increments because they could reflect in a 
different way the presence of coherent structures. 
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FIG. 2. One dimensional longitudinal cut of the velocity 
field v y (xo,y) (continuous line); the same after randomiza- 
tion of the phases (dashed line). In the inset the averaged 
compensated spectrum, with the best fit in the inertial range 
a = 3.26 ± 0.06 

The rationale for this test is to investigate the im- 
portance of correlations between fluctuations at different 
wave-numbers and therefore the "information" content 
brought by coherent structures in 2d turbulent flows. In 
Fig. 2, a one-dimensional cut of the 2-d velocity field be- 
fore and after phases randomization is plotted. At a first 
glance, it is rather difficult to distinguish between the 
true dynamical and the randomized field. This is due to 
the steepness of the spectrum, i.e. only few modes dom- 
inate the real-space configuration. Despite the apparent 
similarity, there are big statistical differences between the 
two fields. Looking at inverse moments we measure a 
clear departure of the true turbulent statistics from the 
bi- fractal prediction (^|), while the randomized configura- 
tions are in agreement with it, with h = 1.13. Due to the 
numerical limited resolution, in order to perform a quan- 
titative statement, one can evaluate only relative scaling 
properties. Therefore, we measure scaling laws of the in- 
verse statistics by plotting all moments T&'iSv) versus 
a reference one, say T( 2 \Sv). This is the same technique 
called ESS |l4|] fruitfully applied in the direct analysis 
of 3d turbulent data with the aim of re-adsorbing some 
finite size effects and extracting scaling information also 



at moderate resolution. Plotting data in this way allows 
quantitative statements only for relative exit-distance ex- 
ponents, x(p)/x(2). 
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FIG. 3. Logarithmic local slopes for the relative scaling 
T^ A \&v) vs T^ 2 \Sv) : for the longitudinal (□) and transverse 
exit-distance (x) for run (A). The same quantities after the 
phases randomization are represented by (A) and (o), respec- 
tively. 

In Fig. 3 we show the logarithmic local slopes for the 
relative scaling of T^(Sv) vs T^(Sv) for longitudinal 
and transversal velocity increments. We stress two im- 
portant results. First, the inverse-statistics moment ex- 
ponents measured on the fields from the numerical sim- 
ulations after phases randomization coincide with the 
prediction (|J) with h = 1.13 as extracted from the av- 
eraged spectrum. On the other hand, the longitudi- 
nal and transversal inverse-statistics moments without 
phases randomization have a more complex, intermittent, 
distribution i.e. they are not described by the bi-fractal 
prediction (^). Second, longitudinal and transversal mo- 
ments are slightly different, indicating that longitudinal 
and transversal velocity fluctuations probe differently the 
smooth part of the 2d field. This is of course connected 
to the fact that longitudinal or transversal velocity dif- 
ferences have different profiles when measured across co- 
herent vortical structures. We also note that transversal 
exit moments display a better scaling behavior than the 
longitudinal ones. 

In Fig. 4, we summarize our results showing the curve 
x(p)/x(2) for both randomized and not-randomized lon- 
gitudinal and transversal exit moments for run (A) fll5| . 
Notice that for p < 1 randomized and not randomized 
data almost coincide because those moments are domi- 
nated by the differentiable fluctuations v(x+r) — v(x) ~ r 
and therefore the relative scaling exponents differ only 
due to the factor %(2) which is almost the same for both 
data sets. On the other hand, clear different statistics 
are measured, for p > 1, by comparing the true dynam- 
ical longitudinal and transversal exit moments with the 
randomized ones. 
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In conclusion, we have given an estimate of inverse- 
statistics moments for signals with a more than smooth 
spectrum, i.e. signals which are differentiable and with 
non-trivial stochastic sub-leading fluctuations. We have 
also shown that statistical properties of a 2d turbulent 
flow are not simply summarized by the spectrum slope. 
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FIG. 4. x(p)/x(2) for 1.5 < p < 6 for both longitudinal 
(A) and transversal (o) exit moments obtained in run (A), 
(□) and (x) indicate the same after randomization, respec- 
tively. Inset: the same but for p < 1. The continuous line is 
the bi-fractal prediction Xbf (p)/Xbf (2), with h — 1.13. Error 
bars have been estimated as the maximum deviation of the 
local slopes from the average value. Run (A) gives indistin- 
guishable results (not shown). 

From the exit-distance analysis it is possible to high- 
light a whole spectrum of more than smooth fluctuations. 
Such fluctuations, being connected with laminar events, 
are the strongest statistical signature of large scale co- 
herence. Experiments with different methods of remov- 
ing/pumping energy at large scales should be performed, 
to investigate the importance of large-scales structures 
in the inverse statistics of flows with different spectra. A 
more quantitative comprehension of the multifractality 
of inverse statistics could be inferred trying to connect 
scaling exponents to the finite time (Lagrangian) Lya- 
punov exponents and the drag coefficient, by extending 
the analysis proposed by Ott and collaborators As 
a final remark we would like to stress that inverse statis- 
tics provide a completely new statistical indicator with 
respect to the standard direct statistics observable. We 
have shown that such method is necessary in all those 
cases where non-trivial fluctuations are sub-leading with 
respect to the differentiable contributions. Obviously, the 
same kind of analysis here reported can be extended to 
temporal signals, opening the possibility of applying the 
method to a broad class of natural phenomena. As an 
example, we just mention possible applications in situ- 
ations common to climatology or meteorology where es- 
timating the probability of persistent velocity configu- 
ration, or of any other dynamical variable, is relevant. 
As a perspective, one important generalization would 



the investigation of multi-dimensional signals by study- 
ing the statistics of d-dimensional volumes between eq- 
uispaced iso-surfaces. The latter method may be, for ex- 
ample, important for analyzing coherence properties of 
two-dimensional or multi-dimensional patterns. 
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